Resonant inelastic x-ray scattering (RIXS) spectra of magnesium diboride 
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Using the tight-binding hnear mufiin-tin orbitals method, the soft x-ray fluorescence K-emission 
spectra of boron in MgB2, excited close to the absorption edge, are estimated. In the calculations 
the angle of incidence 9 between the direction of the incoming photon and the hexagonal axis of the 
specimen is 60° and 75°. Comparison with experiment is possible in the former case where good 
agreement is found. Furthermore, another resonant feature below the Fermi energy is predicted for 
the larger angle. This feature can be related to the excitations to the antibonding B 7r-band in the 
neighbourhood of the L-H line in the Brillouin zone. 

PACS numbers: 71.20-b, 74.25.Jb, 78.70.Ck, 78.70.En 



Magnesium diboride is a new kind of superconductor 
having a critical temperature of 39 K 1]. Although it was 
earlier widely used in chemical technology, the supercon- 
ducting property of this compound was not discovered 
until now At the moment, MgB2 has attracted con- 
siderable attention regarding spectroscopic investigations 
0. For instance, angle resolved photoemission from sin- 
gle crystals have revealed clear dispersions of the occu- 
pied valence bands Additional information has been 
obtained using x-ray absorption - (XAS) and emission 
spectroscopy (XES) |^. In x-ray K-absorption a core- 
hole is created for the Is level by a photon and the re- 
sulting photoelectron is escaping from the specimen. The 
K-emission band can be measured when the deep lying 
level will be filled by a valence-electron and another pho- 
ton is emitted. The resonant inelastic scattering (RIXS) 
differs from the conventional XES in the sense that the 
photoelectron remains in the conduction band above the 
Fermi surface and it will be absorbed by the sample 0. 

Utilizing synchrotron radiation, resonant inelastic x- 
ray scattering is a promising method to probe element- 
specific, local momentum-resolved electronic structure of 
systems that are difficult to investigate using other tech- 
niques. The symmetry of the occupied and unoccupied 
states is coupled with the polarization and direction of 
the incoming and outgoing radiation. Crystal momen- 
tum conservation has been observed in many materials 
e.g. in hexagonal boron nitride and graphite Q. Theoret- 
ical consideration of the RIXS within the band structure 
picture is given e.g. in j^. 

In RIXS, a photon (HiOq) comes to the sample, excites 
a core level (co) electron to the conduction state (cc) and 
a valence electron (e^,) drops to the core hole emitting 
a photon {fuvq'). The formula for the doubly differen- 
tial cross section in the dipole approximation (the wave 
vectors of the ingoing and outgoing photons are small, 
q w g' w 0, compared to the dimensions of the Bril- 
louin zone of the Bloch states {k ) of the electrons) im- 
plemented in our x-ray spectrum program (XSPEC) [oj 
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rris , TTij and t are indexes for conduction 
and valence bands, spin, magnetic quantum number and 
an atom in the unit cell. Here t summation includes 
those atoms in the unit cell which have the core state of 
the specified energy Eq. Using the conventional LMTO 
notation ^lOj the matrix element has the form (IL] 
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where I and m are angular momentum and magnetic 
quantum numbers of the valence and conduction states. 
M^i, Mil ^'^'^ ^rrij I m radial matrix element, its 

energy derivative and angular matrix element, respec- 
tively, having the following properties. 
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where Rt (r) is the radial part of the core-level wave func- 
tion and (l>t,i{r) is that of the valence or conduction band, 
e is the polarization vector of the photon, j refers to the 
total angular momentum of the core state and < j, 
and 1^, m, mg > are the angular parts of the core and va- 
lence/conduction states, respectively. The angular ma- 
trix element leads to the selection rules between the ini- 
tial and final electronic states in the absorption and emis- 
sion of a photon in the scattering process (AZ = ±1). The 
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localized core electron state of an atom t centered at a 
point is 



^ImAr) = Rti\r~rl\)\j,m, > 



(6) 



A'l'l ^ and B^'f ^ are coefficients which depend on the 
crystal structure and the potential of the investigated 
system. The wave function of the valence and conduction 
electrons has the form 



E 
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where 4>ti (r) and ipti (r) are the partial wave and its energy 
derivative. LMTO wave functions are linearized with re- 
spect to energy at some suitable fixed energy E^i. Be- 
cause we are interested in energy eigenvalues E^^j^ near to 

Ei,i we can substitute B^'] coefficients by their approx- 
imate form AW mi^bk ~ ^1^0 J which reduces the required 
memory in computations about a factor of two. 

The electronic structure calculations were performed 
using the scalar-relativistic tight-binding linear muffin- 
tin orbital method in the atomic sphere approximation 
(ASA) . The valence states consisted of Mg 3s, 3p and 
3d states and B 2s, 2p and 3d states. The hexagonal unit 
cell contained three atoms and three empty spheres. We 
used lattice parameters a — b = 2.99 A and c = 3.41 A 
corresponding to equilibrium volume in our calculations. 
For the exchange-correlation potential the parametrized 
form by Perdew and Zunger jl3| was used. The number 
of k points was 648 in the whole Brillouin zone. 

The geometry used in the calculations is the following. 
The direction of the absorbed photon makes an angle 9 
with respect to the hexagonal c-axis of the MgB2 single 
crystal. The absorbed photon is linearly polarized and 
the polarization is in the plane containing the direction 
of the photon and the c-axis of the crystal. In this way 
the component of the polarization vector of the absorbed 
photon along the c-axis can be varied from to 1 by in- 
creasing 6 from 0° to 90°. The direction of the emitted 
photon is perpendicular to the direction of the absorbed 
photon and makes a 90° — 9 angle with respect to the 
c-axis. For the emitted photon we have used two differ- 
ent linear polarizations. The polarization is either in the 
(a, 6)-plane or along the direction of the absorbed pho- 
ton. All our results are expressed using an energy scale 
where the Fermi energy is eV. 

Overall features of both XES and XAS of MgBa are 
well produced by band structure calculations 14]. Thus 
RIXS experiments are also expected to be explained 
within the ordinary band structure picture. Because the 
valence hole is delocalized and the Is core hole in boron 
is well screened jl^ it is reasonable to use ground state 
orbitals to obtain theoretical spectra. Due to the conser- 
vation of both energy and momentum we need detailed 
band structure data in order to interpret the spectra. 
The calculated band structure is shown in Figure^ Since 
we are interested in absorption and emission of photons 
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FIG. 1: Band structure of MgB2. Bands 3, 4 and 5, which 
cross the Fermi level, are shown by a thicker line; solid, 
dashed and dotted, respectively. 
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FIG. 2; RIXS spectra calculated at e = 60°. Excitation 
energy is 0.2 eV relative to the Fermi energy. Polarization 
of the emitted photon is in the (a, 6)-plane (upper spectrum) 
and in the direction of the absorbed photon (lower spectrum) . 
Dashed and dotted lines correspond to spectra calculated us- 
ing bands (3,5) and (4,5), respectively. 



coupled with electronic states near to the Fermi energy, 
the most relevant bands are those which cross the Fermi 
level, i.e. bands number 3-5. Zhang et al. have mea- 
sured RIXS spectra at different angles of incidence {9 — 
15°, 45°, 60° with respect to the c-axis of the MgB2 crys- 
tal) and using different excitation energies (187.25 eV - 
188.25 eV, Fermi energy: Ep = 187.28 eV). They found 
two peaks in their spectra. The intensity of the peak just 
below the Fermi energy depends clearly on both direction 
and energy of the exciting radiation. On the other hand, 
the intensity of the peak at about 2 eV below the Fermi 
energy was considered to be unchanging. 

In Figure |21 we show our results for the 9 = 60° in- 
cidence in the case of two polarizations of the emitted 
photon. The unpolarized emission spectrum is an aver- 
age of these two extreme cases. The peak just below the 
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FIG. 3: RIXS spectra calculated at 61 = 75°. Exitation 
energy (in eV), relative to the Fermi energy, is shown beside 
each spectrum. Polarization of the emitted photon is in the 
(o, fe)-plane. 



ing bands 4 and 5. In the same way, the peak around 
—3.5 eV is mainly due to transitions involving bands 3 
and 5. From FigureQwe see that both of these peaks are 
produced close to the L-H line in the Brillouin zone. In 
this region, near the Fermi energy, the constant energy 
surfaces of the antibonding tt band have a tubular form 

m 

In summary, we have calculated RIXS spectra of 
MgB2 using Kramers-Heisenberg formula with TB- 
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Fermi energy is solely due to absorption to the band 4 
and emission from the band 3. The features between —4 
eV and —1.5 eV consist almost entirely from RIXS in- 
volving band pairs (3,5) and (4,5). By increasing angle 6 
one can measure more directly excitations to the n band. 

In Figure 13 we show the spectra corresponding to 
9 = 75° and the polarization vector of the emitted ra- 
diation, in our calculation, being in the a6-plane. Near 
the top of the valence band there appears a peak which 
disappears when the excitation reaches 0.4 eV above the 
Fermi energy. The broad structure from —4 eV to —1.5 
eV shows not much dispersion up to excitation energies 
0.4 eV above the Fermi energy. However, using excitation 
energies above 1 eV the intensity of the peaks increases 
considerably and the position of the peaks shifts about 
1 eV closer to the Fermi level as the excitation energy 
increases. To decide what bands and what parts of the 
Brillouin zone are responsible for each peak in the spec- 
trum we have calculated some of the spectra band by 
band. As Figure 0| shows the peak around —2 eV (low 
excitation energies) is mainly due to transitions involv- 



FIG. 4: Some of the spectra shown in Fig. 3 calculated using 
the bands (3,4) (dotted), (3,5) (dashed) and (4,5) (solid). 

LMTO eigenvalues and eigenstates for occupied and un- 
occupied electronic states. The results presented show 
that MgB2 can be analysed and interpreted on the basis 
of ground state band structure calculations. In addition 
we predict nonlinear emission structures related with the 
excitation of Is electrons to the unoccupied antibonding 
IT band in the neighbourhood of the L-H line in the Bril- 
louin zone. 
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